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Abstract. We report on several new basic properties of a parabolic dot in the 
presence of a magnetic field. The ratio between the potential strength and the 
Landau level (LL) energy spacing serves as the coupling constant of this problem. 
In the weak coupling limit the energy spectrum in each Hilbert subspace of an 
angular momentum consists of discrete LLs of graphene. In the intermediate coupling 
regime non-resonant states form a closely spaced energy spectrum. We find, counter- 
intuitively, that resonant quasi-boundstatcs of both positive and negative energies exist 
in the spectrum. The presence of resonant quasi-boundstates of negative energies is a 
unique property of massless Dirac fermions. As the strong coupling limit is approached 
resonant and non-resonant states transform into anomalous states, whose probability 
densities develop a narrow peak inside the well and another broad peak under the 
potential barrier. These properties may investigated experimentally by measuring 
optical transition energies that can be described by a scaling function of the coupling 
constant. 



1. Introduction 

Two-dimensional parabolic quantum dots of semiconductor heterostructures have been 
studied widely [TJ E] both experimentally and theoretically because they are excellent 
candidates for single electron transistors. They can effectively confine electrons and the 
number of electrons in them can be controlled using a gate potential. Their Hamiltonian 
is 

H = —(p+ -A) 2 + -mtfr 2 (1) 
2m c 2 

— * — * 

with a magnetic field B is applied perpendicular to the 2D plane (vector potential A is 
given in a symmetric gauge). The characteristic length scale of the problem is given by 
A 2 = — 7 = , where u c = — is the cyclotron frequency. This problem can be solved 

exactly [S], and the eigenenergies are all positive and their spectrum is discrete. 

Massless Dirac electrons [1] moving in a 2D parabolic potential display several 
different features in comparison to massful electrons. They are described by the Dirac 
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Hamiltonian 

e - 1 

H = vpo ■ (p + -A) + -nr 2 . (2) 

C — 

No exact solutions of this problem are known in graphene and several fundamental 
properties are still unknown, such as the existence of resonant and non-resonant states. 
These basic properties may affect the experimentally relevant optical spectrum in a 
profound way. The dimensionless coupling constant of this problem is the ratio between 
the strength of the potential \k1 2 and the LL energy separation Ec = -^r 

a = K f/ Ec = (3) 
nv F 

where the magnetic length is £ = \Jh~c/ Be. One is in the strong coupling regime 
a ^> 1 for small value of B or large value of potential strength k. Parabolic dots in 
magnetic fields have been investigated numerically in the weak coupling regime a < 1. 
The energy spectrum is found to be discrete, and eigenstates are quasi-boundstates 
with long oscillating tails under the barrier [5j El |7J. Also some of these states exhibit 
anticrossings[5j [7J. The parabolic potential acts as a singular perturbation [8] because 
eigenstate wavefunctions are qualitatively different from those in the absence of a 
parabolic potential. 

The problem has not been investigated away from the weak coupling regime. It is 
a highly non-trivial problem. This can be seen as follows. One of the special features 
of graphene LLs is the presence of negative energy states under the potential barrier [7]. 
The first order energy correction of a LL state ip n ,m{r) is, for sufficiently large n, 



(^n,m\V{r)\^ nim ) ~ /t(r 2 ) ~ k1 2 \u\. (4) 

This result suggests that a LL state with a large negative energy, —E C \f2\n\ with 
|n| ^> 1, corresponding to having a large average radius v/ (r 2 ), acquires a significant 
positive energy correction, which can make the renormalized energy positive. In the 
dimensionless units this energy correction is nl 2 \n\/Ec = a\n\, which suggests that even 
for small value of a the correction can be significant for \n\ ^> 1. Moreover, it is unclear 
how eigenstates evolve from weak to strong coupling regimes. A simple dimensional 
analysis suggests that the energy scale of the problem in the strong coupling limit of 
B — > or a — > oo is K 1 / 3 (fryp) 2/ ' 3 - In units of Ec this energy scale is a 1//3 . It indicates 
that the dimensionless energy level spacing increases from ~ 1 to ~ a 1 / 3 as one moves 
from weak to strong coupling regimes. However, studies in ordinary semiconductors 
suggest that the crossover regime may be non-trivial [9]. 

We have investigated these issues by solving large Hamiltonian matrices. Let us 
give a brief summary of our results in a Hilbert subspace of angular momentum J. 
We have studied how eigenvalues and eigenstates evolve as a increase and find that 
they change in a non-trivial way. In the weak coupling limit of a — > the spectrum 
consists of discrete LLs. In the intermediate coupling regime a ~ 1 non-resonant states 
form a closely spaced energy spectrum, see FigJU We find, counter-intuitively, that 
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Figure 1. Schematic energy spectum of a parabolic dot in the intermediate coupling 
regime a ~ 1. Energy spectrum is closely spaced, and resonant states of positive and 
negative energies are present. 



resonant quasi-boundstates of both positive and negative energies exist. The presence 
of resonant quasi-boundstates of negative energies is a unique property of massless 
Dirac fermions, but they are well-defined only for a < 1. In the strong coupling regime 
a ^> 1 both resonant and non-resonant states transform into anomalous states, and a 
sharp distinction between resonant and non-resonant states no longer exists. Probability 
densities of anomalous states develop a narrow peak inside the well and decays slowly 
with small oscillations under the barrier. The energy level spacing between them is 
proportional to the value k}' z [Uvf) 2 ' 3 and is independent of t. We show that optical 
transition energies between resonant quasi-boundstates can be described by a scaling 
function of a. 

2. Basis states and Hamiltonian matrix 

In our Hamiltonian matrix approach the basis states are chosen as graphene LL states 
4>n,m(r) with two components A and B 

WO = C ( - S « n <f*H-'^ ) . (5) 

\ <P\n\,m{r) J 

Here sgn(n) = —1,0,1 for n < 0,n = 0,n > with n and m integers (m > 0), and 
c n = 1 for n = and 1 / \/2 otherwise. These basis states can have positive or negative 
LL energies: 



E n = sgn(n)E cv / 2H- (6) 

The wavefunctions 0„ )m (r*) are the Landau level wavefunctions of ordinary two- 
dimensional systems [TO] 

\ / r\ \m—n\ 



^n,m{r) = A„ im exp yi(n - m)9 - ^ j (J^j 

- \m— n\ ( T 

J {n+m-\m-n\)/2 I 2£2 



U 2 J \£< 

X L\n+m-\m-n\)/2 \ WF?. I ' (^) 
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where A n ^ m are the normalization constants and L™(x) are Laguerre polynomials. In 
the presence of a parabolic potential J = \n\ — m — | remains a good quantum number. 
The average radius of 4> n ,mif) is given by 

(r 2 ) = 2£ 2 (n + m + 1) = 2£ 2 (n + \n\ - J + 1/2). (8) 

To investigate the strong coupling effects a large number of basis states tp n ,m( r ) 
is required. It is convenient to divide the Hilbert space into subspaces of angular 
momentum J = ±|,±|,±|,---. We diagonalize the Hamiltonian matrix in each Hilbert 
subspace J. For given J, the matrix elements of the parabolic potential can be written 
as sum of two components: 



V(r) 

{lpn,m\-E^\^ 



>) 



E c 

c n c n >sgn(nn')(<j)\ 

n\ — l,m 

V(r 



~r C n C 



\n ,m | 



e, 



c 



V(r) 

n'[,m'} • 



101 



n'\ — l,m' 



Using the following property of Laguerre polynomials 



mix) 



- [(n + a + l)L a n -\x) - (n + 1)^1 (*)] 



and the orthogonality 



x a e- x mx)L a (x) dx 



Tin + a + 1) 



(9) 



(10) 



'11' 



TV. 



we evaluate the matrix elements. The resulting matrix is a sparse matrix, see FigfSJ 
The dimension of the Hamiltonian matrix is denoted by N c . When the value of N c is 
sufficiently large the states investigated in this paper do not exhibit dependence on iV c . 



-4-3-2-101234 




Figure 2. Structure of the Hamiltonian matrix is shown for J < 0. For a given J we 
choose to use n as the basis index instead of m. It runs from — (N c — 1) /2 to (N c — 1)/2. 
Filled (empty) squares represent non-zero (zero) elements. For J > allowed basis 
states are given in TABLE II of Ref.[TT]. 
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3. Eigenstates of a Hilbert subspace 

Eigenstates of a Hilbert subspace are obtained by diagonalizing the Hamiltonian matrix. 
They may be written as a linear combination of LL wavefunctions with same angular 
momentum: 



Here quantum number N is chosen to be the value n for which \C n \ is maximum\12\. 

The following exact results [H] are useful in checking numerical results. The value 
of wavefunctions at r = is non-zero only for J = — | and ~: 

WMUj " j finite for J = ±1/2. [U) 

The B and A components of l^^ 1 ^ 2 } and |^/)/ 2 ) are non-zero and are of s-wave type. 
Their values at r = can be written as 



*}&(0) = -^ s g n ( n ) A M-i,M-iCn/V2- (14) 

n 

We will only concentrate on optical transitions involving states with J = 1/2, —1/2, or 
—3/2 since they give the strongest optical strengths. Unless stated otherwise the results 
reported in this paper are for the potential strength k = O.lmeV/nm 2 . 



3.1. Resonant quasi-boundstates 

The computed energy spectrum of J = —1/2 is shown in FigJHJ^a) for the range 
0.029 < a < 0.47 (3.14T < B < 20T). Spectra for other values of J = 1/2 and 
—3/2 are similar, as shown in Sec. IV. Lines labeled by N represent the energies of 
resonant quasi-boundstates and other lines represent non-resonant states. As shown in 
Figj3]^b), in the intermediate coupling regime a ~ 1, the energy levels are closely spaced 
due to negative energy LL states whose energies get strongly perturbed upward by the 
parabolic potential, see Eq.Q. A resonant quasi-boundstate anticrosses other states 
and becomes strongly mixed with the adjacent states, see FigJHJ^b). For example, the 
resonant quasi-boundstate ) is strongly mixed [12] at B = 3.14T (a = 0.47), and, 

as shown in Figj3]^b), there are three states that could be identified as |\l/ ). In these 
states with the energies E = 0.025, 0.019, and O.OlOeV the expansion coefficients C n 
of Eq. f|T2|) take the maximum value at n = with the values Cq = 0.622,0.595, and 

— 1/2\ 

0.287, respectively (Since Co is largest for E = 0.025eV this state is labeled as \^/ )). 
Fig|3]^c) displays probability densities of these states. They form together a resonance 
with the approximate resonant energy E = 0.025eV. In the weak coupling regime the 
width of a resonance is small and, in order to resolve it, the energy level spacing must 
be smaller than the width of a resonance, which requires a large matrix dimension. 
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Figure 3. (a) Eigenenergy spectrum of Hilbert subspace of J = — 4. Lines labeled 
by iV represent resonant quasi-boundstates. 2001 x 2001 matrix is used, (b) Enlarged 
energy spectrum for J = — |. Resonant quasi-boundstate 1 ^ 2 ) anticrosses strongly 
at £> = 3.14T (a — 0.47). Four circles represent these coupled states, (c) Probability 
densities of these four states are displayed. They form together a resonance with the 

— 1/2 

approximate resonant energy e (3.14) = 0.025eV. 




Figure 4. An example of a resonant quasi-boundstate of negative energy (indicated 
as a open circle in Figj3ja)). Dotted (dashed) line represents A (B) components of the 
probability density of l^jj) — l^lj^ 2 ) with energy E — — 0.060eV computed at the 
coupling constant B = 6.88T (a = 0.145). Solid line is the total probability density. 
2001 x 2001 matrix is used. In the absence of the potential the wavefunction has a 
peak near r = (see inset), and its energy is neg ative E = -V2E C = -0.093eV. 
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Note also that resonant quasi-boundstates of negative energies exist. This is a 
unique property of massless Dirac fermions. An example is shown in FigJU They are 
well-defined only for sufficiently large B, i.e., only in the weak coupling regime. The 
appearance of a second peak away from r = in the probability density is different from 
the usual behavior of the wavefunction a resonant state. 

We see that as a increases, or, as B decreases, resonant quasi-boundstates disappear 
into the closely spaced energy spectrum. For N = states this happens around B ~ 3T. 
For larger values of N this happens at smaller values of B. We will show in Sec. Ill (C) 
that, as a increases, the peak at r = increases and the state becomes anomalous. 

3.2. Non-resonant states 




Figure 5. An example of a non- resonant state (indicated as a open circle in 
FigEJa)). Dotted (dashed) line represents A (B) components of the probability density 
\<t> J N ) = |*l} /2 ) with energy E = 0.054eV and B = 2.04T (a = 0.9). As r ->■ the 
dashed line approaches a finite value while dotted line goes to zero. 2001 x 2001 matrix 
is used. 
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Figure 6. Average energy level spacing of non-resonant states near E = as a function 
of N c at B = 1.9T (a = 1). The dimension N c varies from 101 to 18001. Matrix sizes 
N c are 101, 201, 401, 801, 2001, 4001, 9001, 11001, 15001, and 18001. 

The energy spectrum of J = —1/2 in FigJHJ^a) also display non-resonant states. 
Probability density of a non-resonant state looks qualitatively different from that of 
the corresponding unrenormalized LL state. The probability density of a non-resonant 
state |\I/_^ 2 (r)| 2 at B — 2.04T (a = 0.9) is shown in Figj5]^a). Its wavefunction has 
a large peak at r = 0, which is different from the usual behavior of the wavefunction 
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a non-resonant state. As a increases the peak at r = increases even more and the 
state becomes anomalous. Non-resonant states are unique to graphene parabolic wells 
and do not exist in ordinary parabolic wells. In the absence of the parabolic potential 
its energy is E = -V2E C = -0.051eV while in the presence of the potential it is 
e_i (2.04) = l.hEc = 0.054eV. Figj6] displays the energy level spacing of non-resonant 
states as a function of N c at a = 1 (B = 1.9T). We observe that the level spacing 
decreases rather slowly for large N c . However, the energies of resonant quasi-boundstate 
converge rather quickly, see FigJHJ 



3.3. Anomalous states 
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Figure 7. (a) Total probability density of |* 1/2 ) with energy 0.177eV at B = 0.312T 
(a = 15) (solid). Corresponding state at a = is shown as dotted line. The eigenstate 
is obtained by diagonalizing 7001 x 7001 Hamiltonian matrix. Energy levels are also 
shown at B = 0.312T (a = 15) for N c = 7001 and 9001. (b) Expansion coefficients 



C„ = (n|^ 1 ^ 2 ) of anomalous eigenstate |\& 1/z ) 



l/2\ 



The energy spectrum of J = —1/2 in FigJHJ^a) also display anomalous states at low 
magnetic fields of the strong coupling regime a ^> 1. In this regime both resonant and 
non-resonant states transform into anomalous states, and a sharp distinction between 
resonant and non-resonant states no longer exists. Such a state is shown in FigJT^a) 
with the energy E = 0.177eV (in unit of Eq it is 12.6). We see in FigJT^a) that the 
peak value of probability density at r = is much larger than the unperturbed value 
of t^t — 0.16. [13J. For this state the penetration into the barrier should start from 
the turning point r& satisfying \^r\ = E. From this we find that the value r^ft is 
1.3, which is rather different from the estimate of about 0.3 obtained from numerical 
result shown in FigJ7](a) (note that the probability density under the barrier oscillates). 
An anomalous state is a strong coupling effect and can only be obtained correctly by 
computing large Hamiltonian matrices. 
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FigJT^b) displays the expansion coefficient C n as a function of n. Note that C n has 
a long oscillating tail for n < 0. This is intimately related to the probability density 
having a long oscillating tail under the barrier. The sum of C n for n < is approximately 
zero while the sum for n > is finite and makes \l/ (0) large (see Eq. ffT4"l) ). Note that 
the probability density under the barrier is somewhat smaller than that of a = 0. 
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Figure 8. Dimensionless average energy level spacing of the energy spectrum of 
Hilbert subspace J = —1/2 in the strong coupling regime. We have used N c — 9001. 



When B — > or a — > oo the natural length and energy scales of the problem are 
f = (^) 1 /3 anc i K i/3(^ F )2/3 ( Note = a -i/3)_ In units of Ec tnis energy sca i e 

is a 1 / 3 , which should be proportional to the dimensionless energy level spacing of the 
Hilbert subspace of J in the strong coupling regime. Our numerical results in the strong 
coupling regime «>1 are indeed consistent with this, see FigJHJ 

4. Scaling of optical transitions 

4-1- Scaling results 
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Figure 9. Data collapse of dimensionless energies of quasi-boundstate e/ (a) for 
different values of N c and k. Circles, triangles, and squares are for k = 0.1meV/nm 2 . 
Inverted tiangles are for n = 0.2meV/nm 2 . 



In the previous section we showed that energies of the resonant quasi-boundstates 
depend on both k and B. Here we will show that their energies, when measured in 
units of Eq in the limit of large N c , follow a scaling function of a single dimensionless 
variable, namely, the dimensionless coupling constant a, see FigJHJ Fig JIUT a) displays 
the dimensionless energies of resonant quasibound states e J N {a) as a function of a for 
J = —1/2 and 1/2. Fig JTOT b) displays similar results for J = —3/2 and —1/2. 



Resonant, non-resonant, and anomalous states of Dirac electrons in a parabolic well in the presence of magm 




0.0 

Cb)« 
02 

LLJ 

UJ 



-2 
0.0 



21-3/2 




2,-1/2 


1,-3/2 




1,-1/2 


0,-3/2 




0,-1/2 




ii-- 1,-3/2 




♦--1,-1/2 





0.5 1.0 

a 



1.5 



Figure 10 

for J 

J = - 



(a) Energy spectra for J 



I plotted together, (b) Energy spectra 

2001 for 



3 

2 ! 



-i and -| plotted together. Matrix sizes iV c are 2001 for J 
and 2000 for J = \, 



1 

2 • 



Since dimensionless energies of resonant quasi-boundstates e J N (a) satisfy a scaling 
function the transition energies E = ej^,(a) — e J N {a) between them also obey a scaling 

(15) 

This scaling relation holds as long as quasi-boundstates are well defined. 



4-2. Optical transition energies and selection rules 

Before we compute strengths and selection rules of optical transitions let us first mention 
some useful results in computing them. First, in the absence of a parabolic potential 
absorption selection rules are e^ 2 {a) — > e^ 1 (a) for N > and e]^ 2 («) — > e^^(a) 
for iV < 0. These selection rules are displayed schematically in FigJTTJ(see also TABLE 
II in Ref.JTTJ). Due to mixing of different LL states by the parabolic potential the 
selection rule AN = ±1 must be relaxed. Second, the energies of the iV = LL 
states are split by the perturbing parabolic potential and they increase as J decreases: 
e 1 ^ 2 (a) < e 3 ^ 2 (o;) < e Q 5 ^ 2 (a;) < This is also true for other LL states. Third, 
the optical transitions with the strongest strength occur between J = 1/2 and —1/2. 
Other transitions, for example, transitions between J = —3/2 and —1/2 are weaker: the 
transition e (0.25) — > e/ (0.25) has the strength 0.37 with transition energy lAAEc 
while the transition 6q 3 ^ 2 (0.25) — > e 1 1 ^ 2 (0.25) has the strength 0.21 with transition 
energy 1.34Ec (the involved energy levels are shown in FigJTOl). These strengths and 
the corresponding transition energies decrease as the value J of initial states of optical 
transitions decreases. 

Let us use these results to understand what type of optical transitions are possible. 
For this purpose we display in FigJTU for a relatively small value of iV c , eigenenergies 
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Figure 11. We assume that photons are polarized along x-axis. (a) Some examples 
of absorption selection rules for LL states in the absence of a parabolic potential given 
by Eq.©. In fact there are also transitions with the selection rule N — > |AT| + 1 with 
AJ = 1 (N < -1) and N -»• \N\ - 1 with AJ = -1 (N < -2). However, they have 
high energies and are not considered here, (b) Some examples of emission selection 
rules are shown. Selection rules N -> -|JV|-1 with A J = 1 (N > 1) and N -> -|JV|+1 
with A J = — 1 (N > 2) are also possible. These transitions are not relevant here. 
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Figure 12. To display clearly possible optical transitions we plot together, for small 
N c = 101 and 100, the energy spectra for J = — | (black dots) and \ (white dots). The 
value of the strength of the potential K = 0.1meV/nm 2 corresponds to a = 1.424 and 
a = 0.116 at B — 1.5T and B = 8T, respectively. Quasi-boundstatc energies t J N {a) 
are labeled by (N, J). 



for J = — | and | (Results are qualitatively similar to those of a large value of 
N c ). Absorption transition strength between initial and final states and ^/ is 
| (^f\a x \^i) | 2 . Note that there are possible transitions from resonant quasi-boundstates 
of negative energies, which is a new feature. An example of this transition is shown as 
(f) in FigJT2j Its strength is 0.386. Also transitions between resonant quasi-boundstates 
of positive energies are possible. An example is shown as (e) with the strength 0.461. 
Energy of a resonant quasi-boundstate will split into several values at small B due to 
anticrossing with other levels. This will lead to a splitting of transitions. An example is 
shown as (g) and (h). Some examples of absorption transitions involving non- resonant 
states are also shown. For these transitions, due to mixing of different LL states by the 
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parabolic potential, the selection rule AN = ±1 must be relaxed. Examples of these 
transitions are shown as (a), (b), (c), and (d), in FigJT2l Their absorption strengths are 
0.255 (a), 0.316(b), 0.202 (c), and 0.207 (d). 

The next dominant absorption transitions occur for (J, J') = (—3/2, —1/2). In the 
absence of the parabolic potential the selection rules are (a) — >■ (a) for iV > 
and e^ 2 (a) — > %+i(at) for N < 0. Similar selection rules hold for other possible (J, J') 
with AJ = ±1. 

5. Summary 

In each Hilbert subspace of angular momentum J we have studied how eigenvalues and 
eigenstates of a parabolic dot in a magnetic field evolve as a increases. We have found 
that they change in a non-trivial way. In the weak coupling limit of a — > one recovers 
discrete LL spectrum of graphene. In the intermediate coupling regime non-resonant 
states form a closely spaced energy spectrum, see Figj3|a) (The result is different from 
the case of a cylindrical potential, whose the energy spectra in a magnetic field are 
discrete without quasi-b oundst at es [TT| [HI US] except at B = 0[16j). In addition, we 
find, counter-intuitively, that resonant quasi-boundstates of both positive and negative 
energies exist in the spectrum, see Figj3|c). Closely spaced spectrum is consistent 
with the presence of resonant quasi-boundstates [T7] . The presence of resonant quasi- 
boundstates of negative energies is a unique property of massless Dirac fermions, but 
they are well-defined only for a < 1 . In the strong coupling regime of a 3> 1 all resonant 
and non-resonant states become anomalous states that develop a sharp peak in the well 
and decay slowly with small oscillations under the barrier. The average energy level 
spacing in each Hilbert subspace of J approaches a constant value. However, note that 
the total density of states is the sum of each density of states computed in different 
Hilbert subspaces. When a is too large the coupling between K and K' valleys may 
have to be included [18] . 

Optical transitions from resonant quasi-boundstates of positive energies are 
possible. There are also possible transitions from resonant quasi-boundstates of negative 
energies, which is a new feature. Moreover, we find that transition energies between 
resonant quasi-boundstates follow a scaling as a function of the coupling constant 
a. Absorption transitions involving non-resonant states are also possible, and for 
these transitions the selection rule AiV = ±1 must be relaxed due to mixing of 
different LL states by the parabolic potential. It would be also interesting to observe 
experimentally the splitting of optical transition energies due to anticrossing of resonant 
quasi-boundstates with other states. 
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